The Decay of Multiqudit Entanglement 
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We investigate the decay of entanglement of a generalized A-qudit GHZ state with each qudit 
passing through independently in a quantum noisy channel. By studying the time at which the 
entanglement completely vanishes and the time at which the entanglement becomes arbitrarily 
small, we try to find how the robustness of entanglement is influenced by dimension d and the 
number of particles N. 
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I. INTRODUCTION 

Quantum entanglement, as the most non-classical phe- 
nomenon in quantum mechanics, lies in the central posi- 
tion of quantum information theory and has been identi- 
fied as a key resource in many applications such as quan- 
tum teleportation,quantum key distribution and quan- 
tum computation [lj , see Ref. [12j for a review. For large- 
scale quantum information processing, multiparticle en- 
tanglement is indispensable. Therefore the understand- 
ing of the dynamical property of multiparticle entangle- 
ment in realistic environment is of fundamentally impor- 
tance. 

Due to the interaction with environment, multiparti- 
cle entanglement decays inevitably. In the past several 
years, there are many excellent papers concerning with 
the robustness of multiparticle entanglement under the 
influence of environment H & @, 0, B El H3- It 
was shown that for a iV-qubit GHZ state (throughout 
N is even for simplicity) ^(10)®^ + \1)® N ) with each 
qubit exposed to a depolarizing channel independently, 
the bipartite entanglement of (N — n)\n bipartition dis- 
appears in finite time which is named as entanglement 
sudden death (ESD) in Refs. 11 EE El • The entan- 
glement which first disappears is that corresponding to 
l\(N — 1) bipartition and the time at which this entan- 
glement disappears decreases with N. The entanglement 
which last disappears is that corresponding to N/2\N/2 
bipartition and the time at which this entanglement dis- 
appears grows with TV 0, 0] • Aolita et al. [H proposed 
that the time at which the entanglement becomes arbi- 
trarily small is a better quantity characterizing the ro- 
bustness of entanglement than the ESD time and they 
found that this time is inversely proportional to N for 
TV-qubit generalized GHZ state |* 2 ) = a\0)® N + 0\1)® N . 

It is understandable that the robustness of multipar- 
ticle entanglement not only depends on the number of 
particles but also on the dimensionality of each particle's 
Hilbert space. And for different quantum channels, the 
behavior of the multiparticle entanglement will also be 



different. However, it is necessary to have a fixed exam- 
ple to see how all of these results can be obtained. In 
this article, we will study those problems. We suppose 
each particle is in a d-dimensional Hilbert space and it 
is called generally a qudit. We consider a simple case 
of a generalized AT-qudit GHZ state with each qudit in- 
teracts with a quantum channel independently and find 
how does the robustness of entanglement change with the 
increase of d or N through studying the time evolution 
of the bipartite entanglement of (N — n)\n bipartition. 
We use the depolarizing channel and the phase damping 
channel as the quantum channels. 



II. QUANTUM CHANNELS 

First we would like to introduce two d x d matrices 
extremely useful in constructing our quantum channels. 
They are defined as X\i) = \i + l)(mod d) and Z\i) = 
Lo l \i), where i = 0, ..., d— 1 and lo = exp(27ri/d). One 
can easily see that when d = 2, they are just Pauli-sigma 
x and Pauli-sigma z matrix respectively. Now we will 
construct three operator transformations through X and 
Z: 

d-l 

£i{A) = {^-P}A+^Yj X l Z^AZ^X^\ 



ij=0 



d-l 



E 2 (A) = (l-p)A+^J2 ziAzti > W 

where p 6 [0, 1] and A is an arbitrary operator on d- 
dimensional Hilbert spaces. When A is a density ma- 
trix, we can see £i,i — 1,2 as three quantum chan- 
nels. Through calculation, it can be shown that for 
an input density matrix p, £\{p) = (1 — p)p + §1 and 

hip) = (1 - p)p + pJ2k=o Pkk\k)(k\, where 1 is a d x d 
identity matrix. It's obvious to see that in fact £\ is a 
depolarizing channel and £2 is a phase damping channel. 



III. THE EVOLUTION OF ENTANGLEMENT 
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A generalized A^-qudit GHZ state can be written as 
\^d) = J2i=o a i\i)® N 1 where a, is a complex number 



2 



and E*=o \ a i\ 2 — 1- Here we want to show how does 
the bipartite entanglement corresponding to (N — n)\n 
bipartition of it evolve under the influence of the three 
quantum channels constructed above. In this article we 
adopt negativity as the measure of entanglement 

When each qudit is exposed to a depolarizing chan- 
nel, we can calculate ( <3>£Li (l^rf) (^d\) and de- 
note it as Pl (p) ee y££ E£LoE-pKI 2 (§) i 1 - 



d-l 



N-k 



x V 



(1 - P) N Eto Y*Zli*i <Xia*j(\i){j\r N , where V means 
all possible permutations. Partial transposing the part of 
n particles of p\ (p) and noting that the first term of p\ (p) 
is diagonal so it will not be changed after partial trans- 
pose, we have pi(p) r ™ = J2t=o EfeLo J2v H 2 (l) 



1 - 



d-l. 



-P 



N-k 



x V 



\\i)(i\r( N -v®(EU,&i li)01J 
(i -p) N EtoE-Zo^ a ^^ {N ~ n) ® Wl)®"- 

There are d ^ d ~ 1 '> eigenvalues p of pi(p) Tn that can 



be negative and are 
eigenvalues of the following 



determined 

d(d-l) 
2 



by the smaller 
2x2 matrices 



a:/ 



n 



ana* {l -p) 



N 



a*a 3 {l-p) N A|_„ 



where i < j and A*f 

N — 

{ I _ i_ 

d 



1 - 



(\ Ti / \ iV — Tt / \ -i 

5) C 1 - 1 ^) + Kf(s) 

One can easily derive /i^ = — J (£ r Y) 2 — rfn , where 
C = |(A« + A^_„) and r$ = A«Ajj_ - \a t a,\ 2 (l - 
p) 2N . Then we can define J\[$ = max{— 0} and the 



negativity can be obtained by summation J\f n {p\{p)) = 
d-i 

i<j = 



Similar to the case of qubit 0, Aff < Af 2 3 < ... < 
^ 0r a Si yen P & i r °f a i ancl a ji which straightfor- 
wardly leads to M\ < N2 < • •• < A/jv/2- So the bipartite 
entanglement corresponding to the most balanced par- 
tition still disappears last whereas the one correspond- 
ing to the least balanced partition disappears first. Let 
p n = max^j p 1 ^ , where p 1 ^ is the solution of the equation 
pl{ = 0, then N n {pi{pn)) = 0. Now we investigate dy- 
namical property of Af N / 2 , which disappears last. First 
we want to know when does it completely vanish. By 
solving ls N / 2 = 0, it's easy to find that 



2|a{Qij| N d 



'N/2 



{\0Li 



,\ 2 Y- 



y 



(2) 



which coincides with Eq.(6) of Ref. Q when d — 2. After 
p reaches the value pn/2 = max i<j Pjv/2' ■A/iv/2 = 0. It's 
obvious that p N / 2 < 1 so Pn/2 < 1 an d so ESD happens. 
Before p — Pm/ 2 , A/jv/2 > and pi{p) must be still entan- 
gled. Second we want to know when does AFn/2 become 
arbitrarily small, which is practically important because 
before the entanglement is zero, it can be so small that 
it's useless as a resource. Like Rcf48], we suppose e is 
an arbitrarily small positive number and define a criti- 
cal probability p l J such that p l ^ 2 {p^) = ep l ^ 2 (0), which 
leads to an equation 

-\a iaj \(l - p\ 3 ) N = -e\a iaj \. (3) 

When p > p e ee maxi <:) - pl J , we can think A/jv/2 is too 
small to be used as a resource. 

As the second case let's consider the situation 
where each qudit is exposed to a phase damp- 
ing channel. Similar to the case of depolariz- 
ing channel, it can be obtained that p 2 {p) = 

(0ti^)(I^X^I) = EtoKI 2 (K)M)^ + (i- 

P) N Ef=o Ej=oj^ aia*{\i)(i\)® N and the partial trans- 



pose of it is p 2 {pf- = T,to\ai\ 2 (\i)(i\)® N + (1 - 

rf^EtoES^i^iCWO'l) 8 ^ ® (\Mi\r n - so 

it's easy to verify that the negativity M n {p 2 [p)) can 
be expressed as E»</=o max {0i — v n '}> where f^- 7 = 
— \aiaj\(l—p) N is independent of n. We note that for any 
n, N n (p 2 (p)) — only when p = 1, meaning that no ESD 
happens for phase damping channel. When the negativ- 
ity becomes arbitrarily small, namely v^[j> % [) = ei^ (0), 
we have an equation (1 — pf) N = e. 

In what follows we will focus our attention on the re- 
sults of entanglement evolution under the depolarizing 
channel and phase damping channel to study how does 
the robustness of entanglement change with N and d. 



IV. ROBUSTNESS OF ENTANGLEMENT 

First we fix d to study the relation between the en- 
tanglement robustness and N. For the depolarizing 
channel, Eq.@ tells us that the ESD time of Njy/2 
grows with N (FIG.l). Noting that YrniN-^ooP^^ = 
2d/ (2d+ 1 + \/5) which is independent of i and j, we find 
Pn/2 — maxj<j p % ^ 2 becomes closer to this value while N 
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FIG. 1: (color online) For the depolarizing channel, we 
demonstrate how can pjv/2 be influenced by both N and d. 
Here we choose a; = 1/Vo for convenience and N = 4 (black 
cubic), N = 6 (red circle) and N — 8 (blue triangle) respec- 
tively, (a): the behavior of Pjv/2 for large d. (b): the behavior 
°f Pjv/2 for small d. It's easy to see that Pn/2 grows with N 
for a fixed value of d saturating to 2d/(2d + 1 + VE) and it 
also grows with d and saturates to 1 for any fixed N. 



FIG. 2: (color online) For the depolarizing channel, we show 
the relation of p e with N and d. Here we choose oti = 1/ \/d 
for convenience, e = 0.01 and N = 4 (black cubic), AT = 6 
(red circle) and N — 8 (blue triangle) respectively, (a): the 
behavior of p £ for large d. (b): the behavior of p t for small 
d. It can be seen that when d is fixed, p e decreases with TV 
while when N is fixed, it grows with d to a saturated value 



grows (FIG 1). Moreover, when N is big enough, consid- 
ering p 1 / is very small, the first term in the LHS of Eq. ((3]) 
can be omitted, leading to p t ~ —(l/N) lne (FIG. 2). For 
the phase damping channel, no matter what N is, no 
ESD happens. However, p e ~ —(l/N) lne still holds. We 
can see that so long as d is fixed, the scaling relation be- 
tween p e and N is always the same with that in Ref . [|[ , 
where d = 2. 

Then we fix N to study the relation between the en- 
tanglement robustness and d. For the depolarizing chan- 
nel, from Eqs.© and one can find that when the 
coefficients a«s are given (it's not necessary to set all 
ajS the same), both p^/ 2 (also Ptq/2) an d p l J (also p e ) 
grow with d, meaning the entanglement is more robust 
(FIG.l and 2). Moreover, when d is large enough, we 
have lim^oo p N / 2 = 1 (FIG.l) and lim^oo p e = l-e 1/N 
(FIG. 2). For the phase damping channel, when N is 
fixed, it's obvious that both the ESD time (infinity) and 
p l J — 1 — e 1 /^ are independent of d. 

One question arises that for the depolarizing channel, 
what's the dynamical behavior of the bipartite entangle- 
ment corresponding to the least balanced partition that 
vanishes first? As we know, in qubit case (d=2) it van- 
ishes earlier when iV grows. Now our numerical calcula- 
tion shows its relation with N and d (FIG. 3). 
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FIG. 3: (color online) For the depolarizing channel, we show 
the relation of pi with d and N. Here we choose cu = 1/ yd 
for convenience and N = 4 (black cubic), N = G (red circle) 
and N — 8 (blue triangle) respectively. It can be seen that 
when d is fixed, pi decreases with TV while when JV is fixed, 
it grows with d. 
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One could expect that the dependence of entanglement 
robustness on N and d would be similar, since in both 
cases we are increasing the dimension of the Hilbert space 
of each of the bi-partitions. But in fact their influences 
on the entanglement robustness are different. Roughly 
speaking, entanglement increases with d and thus it is 
more robust, while entanglement becomes more fragile 
with N since more particle are entangled together and 
it becomes easier to be destroyed. We take the depo- 
larizing channel as an example, for which according to 
the discussion above, p e increases with d while it de- 
creases with N. This can be explained as follows. If 
d is fixed, when we increase N, we increase the compo- 
nents of the state. Considering the depolarizing channel 
acts locally on every component, the growth of N will 
make the entanglement more fragile. In the depolariz- 
ing channel £i(p) — (1 — p)p + §1, p can be regarded as 
the probability with which p is broken by the channel. 
The probability of the iV-particle state to be broken by 
the collective local action of N channels on each particle 
must grow with N (about Np as a rough estimation), 
leading that the entanglement becomes less robust. The 
probability with which the TV-particle state is not broken 
can be estimated roughly as (1 — p) N and this probability 
can also be expressed as e considering the entanglement 
decay. Therefore we have (1 — p) N — e, leading to our 
familiar results. If N is fixed, with the increase of d, 
the influence of the channel on p will be more and more 
negligible, which leads the entanglement becomes robust. 

V. SUMMARY 

In this brief report, we mainly investigate the dynam- 
ical property of entanglement of a generalized ./V-qudit 



GHZ state under the influence of the depolarizing channel 
and the phase damping channel. We study the relation of 
the entanglement robustness with N and d. First we con- 
sider the ESD time t\ and respectively of the bipartite 
entanglement corresponding to the most balanced parti- 
tion (t\) and the bipartite entanglement corresponding 
to the least balanced partition [t-z). When d is fixed, for 
the depolarizing channel t\ delays when N grows whereas 
t2 becomes earlier. For the phase damping channel, no 
ESD happens for any N. These results are qualitatively 
the same with that of d = 2. When N is fixed, for the 
depolarizing channel both t\ and ti grow with d whereas 
still no ESD happens for the phase damping channel for 
any d. Next we consider the time at which the bipar- 
tite entanglement becomes arbitrarily small (£3). When 
d is fixed the scaling relation between p e and N is to- 
tally independent of d therefore the same with the result 
in Ref.[|[ for both channels. When N is fixed, for the 
depolarizing channel £3 grows with d whereas it's inde- 
pendent of d for the phase damping channel. There are 
many other multiqudit entangled states such as general- 
ized 7V-qudit W-state [9J and other quantum channels. It 
is not clear how the behaviors of the entanglement differ 
in those situations. This is worth being studied further. 
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